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1.  INTRODUCTION 

On  7  February  1983,  the  Office  of  Naval  Research  (ONR)  convened  a  Workshop 
on  Low-Reynolds-Number  Aerodynamics  (in  the  Reynolds  number  range  from  200,000 
to  1,000,000)  to  review  and  appraise  the  status  of  aerodynamic  data  and  theory 
for  the  design  of  low-speed,  low-altitude  remotely  piloted  vehicles  (RPVs).^ 
Such  RPVs  are  potentially  to  be  deployed  as  decoys  of  a  naval  fleet,  confusing 
the  enemy's  radar  system  so  that  the  positions  of  the  fleet  cannot  be  accu¬ 
rately  detected.  A  special  focus  program  was  subsequently  established  at  ONR 
to  investigate  low-Reynolds-number  and  low-altitude  aerodynamics.  The  present 
report  presents  our  effort  in  investigating  the  effects  of  atmospheric  un¬ 
steadiness  on  the  performance  of  a  full-scale  Wortmann  FX-63-137  wing. 

The  status  of  the  knowledge  of  aerodynamics  for  low-speed,  low-altitude 
decoys  can  be  reviewed  conveniently  from  several  technical  viewpoints:  an  un¬ 
perturbed  versus  a  real  atmosphere;  two-  versus  three-dimensional  flows;  and 
lift  coefficients  less  than  approximately  1.5  to  1.7  versus  higher  lift  coeffi¬ 
cients.  To  date,  essentially  all  the  research  data  have  been  obtained  in  what 
was  hoped  to  be  an  unperturbed  flow.  Unfortunately,  the  data  are  contaminated 

to  some  degree  by  disturbances  inherent  in  the  test  procedure  (wind  tunnel  or 

1  2 

flight).  Based  on  the  most  reliable  data,  *  at  least  for  the  higher  Reynolds 
number  flows,  lift  coefficients  up  to  roughly  1.5  to  1.7  are  obtainable  for 
two-dimensional  airfoils  ("unperturbed  atmosphere”)  designed  for  high  lift 
that  have  been  experimentally  verified  for  a  limited  number  of  cases. 

3 

For  the  lower  range  of  Reynolds  numbers  of  interest,  to  quote  Lissaman, 

”.  .  .the  design  of  low-Reynolds-number  airfoils,  particularly  at  Reynolds 

a 

numbers  below  300,000,  is  still  a  black  art"  (also  see  Mueller  ).  A  qualita¬ 
tive  picture  of  the  flow  around  an  airfoil  with  lift  coefficients  in  this  range, 
under  steady  conditions,  is  available.  Detailed  quantitative  information  on 
the  upper  surface  bubble  structure,  its  role  in  the  transition  process,  and 
the  transition  process  itself  is  not  available,  especially  for  an  unsteady  en¬ 
vironment  such  as  the  lower  atmospheric  boundary  layer  over  the  ocean  surface. 
Although  there  is  a  dearth  of  detailed  pressure  data  for  the  three-dimensional 
wings,  Lieback  has  commented  that  rather  simplified  design  procedures  have  met 
with  reasonable  success  in  predicting  the  gross  aerodynamic  characteristics  of 
three-dimensional  wings. 

A  lift  coefficient  in  the  range  of  1.5  to  1.7  is  quite  low  for  even  ap¬ 
proaching  a  down-graded  compromise  in  the  RPV  requirements . ^  It  appears  that 
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any  attenpts  to  obtained  significantly  higher  lift  coefficients  will  require 
more  exotic  airfoils  (fixed  multi-element  airfoils,  movable  multi-element 
airfoils^  or  blowing).  A  compromise  among  the  design  criteria  is  clearly 
required  when  such  added  complexity  is  incorporated  into  the  design  of  future 
RPVs.  In  any  case,  the  development  of  such  higher  lift  airfoils  would  benefit 
by  a  better  understanding  of  the  flow — an  understanding  that  can  be  achieved 
only  through  extensive  experimentation.  An  extensive  experimental  program  is 
beyond  the  scope  of  the  present  work.  Instead,  the  present  experiments  are 
dedicated  to  investigating  one  of  the  major  uncertainties  of  the  performance 
of  low-Reynolds-number  wings,  i.e.,  the  effects  of  atmospheric  unsteadiness. 
While  such  effects  are  known  to  be  significant  in  many  instances no  method 
of  quantification  exists.  Even  if  the  atmospheric  turbulence  and  related 
parameters  were  known  in  detail  (though  according  to  Spence1  they  certainly 
are  not  available),  their  impact  on  the  aerodynamic  performance  of  a  wing 
under  such  conditions  could  be  estimated  only  in  the  crudest  fashion. 

For  a  low-speed  and  low-altitude  decoy  RPV,  the  gust  and  turbulence  In  the 
atmospheric  boundary  layer  are  expected  to  have  significant  effects  on  its  per 
formance.  One  of  these  effects  is  due  to  the  unsteadiness  of  the  wind  field, 
which  effectively  causes  the  wing  of  the  RPV  to  experience  unsteady  pitching/ 
plunging  notions  when  the  reference  frame  is  fixed  on  the  wing.  In  particular 
under  the  Influence  of  the  low-frequency  fluctuating  wind  components,  the  wing 
is  predominantly  subject  to  an  unsteady  plunging  rather  than  pitching  motion 
relative  to  the  Instantaneous  wind  field.  When  an  airfoil  undergoes  oscil¬ 
latory  motions,  i.e.,  pitching  or  plunging,  several  dramatic  phenomena  that 
deviate  greatly  form  the  static  situation  take  place.  A  brief  review  of  the 
unsteady  aerodynamics  is  given  below.  For  more  detailed  information  relevant 

g 

to  the  present  investigation,  the  readers  are  referred  to  Liu. 

When  the  angle  of  attack  oscillates  around  a  mean  value  a  near  or  greater 

o 

than  the  static  stall  angle,  the  maximum  values  of  lift,  drag,  and  pitching 
moment  (C^,  Cp,  and  |d^|,  respectively)  can  greatly  exceed  their  static 
counterparts.  The  higher  value  in  the  lift  coefficient  is  generally  referred 
to  as  the  "lift  overshoot."  Such  increases  in  CT ,  C_.,  and  |cw|  cannot  be 

1»  U  M 

reproduced  even  qualitatively  by  neglecting  the  unsteady  notion  of  the  airfoil 
During  the  oscillatory  cycles,  large  hysteresis  develops  in  the  aerodynamic 
forces  and  moments  with  respect  to  the  instantaneous  angle.  The  above  process 
is  referred  to  as  dynamic  stall,  which  is  intimately  related  to  unsteady 


TR-43 1/11-87 


separation  of  the  flow  from  the  oscillating  airfoil.  The  process  of  dynamic 

stall  may  be  classified  into  three  important  regimes,  namely,  stall  onset, 

light  stall  and  deep  stall,  depending  on  the  state  of  the  viscous-inviscid 

9  10 

interaction  for  an  oscillating  airfoil.  *  Stall  onset  is  the  limiting  case 
of  the  maximum  unsteady  lift  that  can  be  obtained  with  no  significant  penalty 
in  pitching  moment  and  drag.  Further  increase  in  a  beyond  results  in 

dynamic  stall,  which  is  characterized  by  large  phase  lags  and  hysteresis  in  the 
separation  and  reattachment  of  the  viscous  flow.  In  the  light  stall  regime, 
the  scale  of  the  flow  separation  remains  on  the  order  of  the  airfoil  thickness. 
As  the  maximum  incidence  increases  to  values  well  beyond  the  static  stall 
angle,  the  flow  over  the  airfoil  is  in  the  deep  stall  regime.  A  strong  vortex¬ 
like  disturbance  initiates  in  the  leading  edge  region  and  subsequently  sheds 
from  the  airfoil  and  moves  downstream  over  the  upper  surface  of  the  airfoil. 

The  process  produces  large-amplitude  hysteresis  in  CT ,  CL,  and  Cw.  The  scale 
of  the  viscous  interaction  zone  is  on  the  order  of  the  airfoil  chord. 

For  an  airfoil  undergoing  a  periodic  plunging  motion,  a  similar  response  in 

the  deep  stall  regime  has  been  shown, ^  but  significantly  different  responses 

12  13 

in  the  light  stall  regime  have  been  reported.  *  In  essence,  the  stall  cell 
Is  not  well  organized  and,  when  it  occurs,  it  has  no  orderly  downstream  propa¬ 
gation.  The  vortex-like  disturbance  moves  downstream  over  the  airfoil  surface 
at  a  relatively  slow  speed.  The  vortex  does  not  have  a  strong  peak,  and  the 
low-pressure  zone  formed  on  the  upper  surface  is  elongated.  As  a  result,  a  high 
lift  coefficient  could  be  sustained  for  a  relatively  long  period  of  tine  when 
the  angle  of  attack  increases  beyond  the  static  stall  angle. ^ 

The  studies  of  Maresca  et  al.^  also  established  the  existence  of  vortex- 
shedding  phenomena  and  strong  hysteresis  effects  due  to  translational  oscilla¬ 
tions  in  the  streamwise  direction  and  at  various  oblique  angles  relative  to  the 

oncoming  flow.  Unsteady  effects  were  found  to  be  strong  enough  to  prevent 

14 

stalling  during  part  of  the  cycle.  St.  Hilaire  et  al.  examined  the  effects 
of  wing  sweep  at  an  angle  X  with  respect  to  the  free  stream  for  an  oscillating 
wing.  Their  results  show  that  the  wing  sweep  tends  to  delay  the  onset  of  dyna¬ 
mic  stall  and  to  reduce  the  rate  of  change  of  and  as  the  stall 
begins.  It  also  reduces  somewhat  the  amplitude  of  the  hysteresis  loops. 

Another  unsteady  phenomenon  of  great  concern  pertains  to  the  aerodynamic 
damping.  For  a  wing  undergoing  unsteady  oscillation,  the  instantaneous  work 
done  on  the  fluid  by  the  wing  due  to  its  motion  is  dW  ■  -Mda,  where  M  is  the 
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pitching  moment  about  the  axis  of  the  rotation  and  is  positive  for  nose-up  air¬ 
loads.  Normally,  C  da  is  negative,  but  it  can  become  positive  during  some 

M  10 

phase  of  dynamic  stall.  As  a  result,  work  is  being  done  on  the  wing  by  the 
fluid.  A  pitching-damping  parameter,  defined  as  ;  -  /  Cwda/4a,  is  a  measure  of 
the  net  work  done  over  a  complete  cycle  of  unsteady  motion.  If  5  is  negative, 
the  wing  extracts  energy  from  the  airstream,  and  the  unsteady  motion  may  in¬ 
crease  in  amplitude  unless  restrained.  This,  of  course,  is  the  condition  of 
flutter.  Stall  flutter,  arising  from  this  negative  pitch  damping,  tends  to 
occur  when  the  airfoil  is  oscillating  in  and  out  of  stall.  The  phenomenon  of 
stall  flutter  is  illustrated  in  Figure  i.^,10  ^  shown  in  the  figure,  x.  is 

positive  throughout  the  cycle  when  a  =7.3°  (not  stalled)  and  when  a  =  24.6° 

"  o  o 

(remain  stalled).  For  aQ  “  14.9°,  however,  the  airfoil  is  stalled  part  of  the 
time  and  not  at  all  at  other  times,  and  the  areas  inside  the  clockwise  loop 
( 5  <  0)  and  counterclockwise  loops  ( 5  >  0)  are  approximately  the  sane,  indi¬ 
cating  neutral  stability.  Had  the  mean  angle  been  slightly  less,  the  net 
damping  would  probably  have  been  negative,  and  the  oscillation  would  therefore 
have  been  unstable  if  unconstrained.  Stall  flutter  may  also  occur  In  plunging 
oscillation  when  the  wing  is  in  and  out  of  stall  under  the  influence  of  a  rela¬ 
tively  large-amplitude  fluctuation  of  the  vertical  wind  components,  which  may 
potentially  exist  near  the  ocean  surface  at  high  sea  states.  Stall  flutter  may 
lead  to  structural  damage  due  to  material  fatigue  provided  the  frequency  of  the 
stall  flutter  Is  in  resonance  with  the  natural  frequency  of  the  RPV  wing. 

To  fill  this  serious  gap  in  our  knowledge  of  the  environmental  effects, 

Flow  Research,  Inc.  (FLOW)  has  engaged  in  a  unique  experimental  investigation 
to  quantify  such  effects  on  the  low-Reynolds-number  aerodynamic  characteristics 
of  full-scale  wings  representative  of  an  advanced  RPV  configuration.  The  main 
thrust  of  FLOW'S  involvement  was  to  collect  and  analyze  full-scale  experimental 
data  and  to  understand  the  behavior  of  the  flow  field  around  low-Reynolds- 
number  (LRN)  wings  under  various  environmental  conditions  often  encountered  by 
RPVs  during  their  operations.  Up  to  the  present,  we  have  investigated  the 
effects  of  atmospheric  gust  and  turbulence.  The  efforts  were  divided  into  two 
phases.  During  Phase  I,  an  environmental  aerodynamic  test  system  (EATS)  was 
designed  and  built. ^  In  Phase  II,  several  series  of  experiments  were  conduc¬ 
ted  to  Investigate  the  effects  of  atmospheric  unsteadiness  on  the  performance 

of  the  Wortraann  wing.^  The  results  of  full-scale  testing  were  compared  with 
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the  wind  tunnel  data  obtained  with  the  same  wing  under  steady  conditions.  ’ 
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The  wind  tunnel  results  serve  as  baseline  cases  for  assessing  the  unsteady 
effects  of  the  wind  field  on  the  performance  of  the  Wortmann  wing.  Much  of 
the  experimental  results  and  their  interpretation  and  assessment  have  been 
reported  elsewhere.® 

In  this  report,  we  summarize  the  most  recent  work  performed  in  the  research 

8  16 

project.  Most  of  the  work  has  been  reported  in  two  technical  papers,  ’  and 
frequent  reference  to  these  papers  is  made  to  avoid  unnecessary  duplication  and 
to  ensure  continuity.  In  addition  to  the  results  reported  in  these  papers, 
the  analyses  of  the  high-frequency  characteristics  of  the  wind  components  and 
force/moment  are  presented  in  this  report  to  explore  other  potentially  impor¬ 
tant  effects  due  to  the  unsteady  nature  of  the  wind  field.  We  present  the 
time  series,  the  hysteresis  and  the  spectra  of  the  ambient  and  aerodynamic 
parameters.  In  particular,  we  attempt  to  explore  the  dynamic  stall  phenomenon 
that  might  lead  to  stall  flutter.  Stall  flutter  must  be  considered  in  the 
design  and  fabrication  of  the  wings  of  RPVs.  If  not  controlled,  this  flutter 
nay  result  in  serious  structural  damage  to  the  RPV,  particularly  when  complex 
wing  shapes,  such  as  a  multi-element  wing,  are  used  to  increase  aerodynamic 
performance  at  the  expense  of  structural  integrity. 
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2.  ENVIRONMENTAL  AERODYNAMIC  TEST  SYSTEM  (EATS) 

During  Phase  I,  an  environmental  aerodynamic  test  system  (EATS)  was  designed 
and  assembled  for  testing  full-scale  wings  under  a  broad  range  of  atmospheric 
conditions,  including  gust  and  turbulence,  cross  winds,  rain  and  possibly  snow. 
In  its  present  form,  the  EATS  consists  of  a  full-scale  Wortmann  FX-63-137  wing 
(smoothed  by  Eppler)  with  a  span  of  3.67  m  and  a  chord  of  0.61  m  mounted  on  an 
instrumented  moving  platform,  a  one-ton  Ford  truck.  The  EATS  is  equipped  with 
a  three-component  aerodynamic  balance  (lift,  drag  and  pitching  moment)  mounted 
inside  the  wing  body  and  on  a  boom  about  6-ra  above  the  ground.  In  addition  to 
the  forces  and  moment,  the  EATS  is  capable  of  simultaneous  measurements  of  the 
atmospheric  parameters ,  such  as  the  three  wind  velocity  components ,  the  baro¬ 
metric  pressure  and  temperature.  For  a  detailed  description  of  the  EATS  and 
its  operation,  refer  to  References  8  and  16. 

The  EATS  was  developed  based  on  the  design  of  a  similar  system  for  testing 
19 

hang  gliders.  For  example,  the  original  design  of  the  flexures  used  in  the 

aerodynamic  balance  were  essentially  a  scaled-down  version  of  those  used  by 
19 

Pregler  et  al.  Miniature  load  cells  with  superior  performance  characteris¬ 
tics  (in  terms  of  sensitivity  and  temperature  drift)  to  those  used  by  Pregler 
et  al.  were  incorporated  into  the  new  balance.  Extensive  tests  were  conducted 
to  determine  the  performance  of  the  EATS.  Based  on  the  test  results,  modifi¬ 
cations  of  a  few  components  were  made  to  improve  the  performance  characteris¬ 
tics.  In  particular,  the  scaled-down  flexures  were  found  to  be  too  weak  to 
handle  the  load  generated  by  the  full-scale  Wortmann  wing.  The  flexures  were 
sheared  off  repeatedly  at  the  narrow  necks  during  testing.  Subsequently,  the 
flexures  were  redesigned  using  beryllium  copper  with  a  thin  straight  section 
2.54  cm  long  and  0.24  cm  in  diameter,  instead  of  using  high-strength  tool  steel 
with  two  short  narrow  necks.  The  improved  flexures  worked  reliably  without  a 
single  failure  during  all  the  experiments.  Another  modification  was  to  add 
two  flexures  (without  the  load  cells)  to  increase  the  stability  in  the  yaw 
direction.  But  this  resulted  in  a  loss  of  the  capability  of  providing  the  yaw 
measurement.  To  match  the  dynamic  response  of  the  balance  and  the  wind  sensors 
(Gill  propeller  anemometers  with  a  dynamic  response  of  2  Hz  at  8  m/s),  4-  to 
5-Hz  low-pass  filters  were  installed  in  the  output  circuitries  of  the  forces 
and  moment  controller. 

Two  mini-clinometers  (Sperry  Corporation)  were  used  to  set  and  monitor  the 


angle  of  attack  and  the  leveling  of  the  truck  in  the  two  horizontal  planes. 


C*' 
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Leveling  of  the  EATS  in  the  horizontal  plane  was  achieved  by  jacKing  up  the 
truck  at  three  points  under  the  main  body  frame,  once  the  prevailing  wind 
direction  was  determined.  A  third  clinometer  was  used  to  monitor  the  inclina¬ 
tion  angle  of  the  Gill  anemometers. 

During  the  initial  tests,  we  discovered  that  the  force/moment  data  were 
seriously  contaminated  by  the  vehicle  vibration  induced  by  road  roughness  when 
the  EATS  was  operating  in  a  run  mode  along  a  runway  or  an  automobile  racing 
track.  Significant  modifications  to  incorporate  a  sophisticated  damping  system 
would  be  required  to  remove  the  contamination.  These  modifications  were  not 
possible  due  to  budgetary  constraints.  The  run-mode  approach  was  then  dis¬ 
carded  for  the  present  investigation.  To  avoid  vehicle-induced  contamination, 
several  series  of  experiments  were  conducted  under  a  park-mode  configuration 
by  directing  the  wing  into  the  prevailing  wind. 

The  results  were  grouped  into  two  series  with  Reynolds  numbers  in  the  range 
of  80,000  to  200,000  and  250,000  to  450,000.  The  first  series  (low  R^)  was 
conducted  at  the  parking  lot  of  FLOW  during  windy  days.  Depending  on  the  wind 
direction,  the  upwind  ground  features  were  either  a  large  vacant  lot  with  low 
brush  (facing  south)  or  a  long  alley  formed  by  several  warehouse  buildings 
(facing  north).  The  wind  field  was  quite  unsteady  with  the  wind  speed  and 
direction  varying  noticeably  during  the  experiments.  When  the  wind  direction 
changed  momentarily,  vortices  shed  from  the  corners  of  buildings  and  tree  tops 
were  observed  from  a  flow  visualization  experiment  using  smoke  as  the  tracer. 

The  reference  runs  were  conducted  downwind  of  a  warehouse  building.  The 
results  were  carefully  screened  to  discard  those  runs  during  which  significant 
changes  in  the  wind  conditions  took  place. 

The  second  series  of  experiments  (high  R£)  was  conducted  at  the  Ellensburg 
Airport,  Ellensburg,  Washington,  located  east  of  Snoqualmie  Pass  in  the  Cascade 
Mountains.  During  the  months  of  March  through  July,  steady  and  strong  seasonal 
winds  up  to  20  m/s  are  common  through  the  pass.  The  upwind  terrain  is  excep¬ 
tionally  flat  for  miles  with  very  low  desert  plants.  During  the  experiments,  the 
wind  conditions  were  relatively  steady  except  for  the  turbulence  generated  by 
the  flat  terrain  upwind.  The  reference  runs  were  conducted  in  a  hanger,  which 
provided  an  ideal  and  quiescent  environment.  For  each  angle  of  attack,  we  con¬ 
ducted  a  series  of  about  10  runs.  Reference  runs  were  conducted  before  and/or 
after  each  series  of  runs.  See  Reference  8  for  a  description  of  the  procedures 
for  setting  up  the  reference  runs  and  for  the  force/moment  measurements. 
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As  reported  by  Liu,  the  quality  of  the  data  obtained  in  the  second  series 
was  significantly  better  than  that  of  the  first  series.  In  particular,  the 
data  scatter  is  noticeably  higher  for  the  low-Reynolds-number  series  than  for 
the  high-Reynolds-number  series,  as  expected. 
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3.  RESULTS  AND  DISCUSSION 

Most  of  the  essential  results  were  reported  previously  in  the  papers  by  Liu 
16  fl 

et  al.  and  Liu.  The  first  paper  presented  a  detailed  description  of  the 
EATS  and  some  preliminary  test  results.  The  second  paper  discussed  the  latest 
modifications  of  the  EATS  and  presented  the  results  of  the  two  series  of 
full-scale  experiments.  Refer  to  these  two  papers  for  detailed  information. 

In  this  report,  we  present  additional  analyses  that  are  oriented  toward 

investigating  the  unsteady  nature  of  the  results.  First,  we  present  four  sets 

of  time  series  covering  a  range  of  mean  angles  of  attack  from  small  to  beyond 

dynamic  stall.  Second,  we  examine  the  unsteady  aerodynamics  by  selecting  from 

the  time  series  several  individual  cycles  with  large  amplitudes.  In  particular, 

we  illustrate  the  unsteady  loops  of  C^,  Cp,  and  for  three  different  cycles 

for  angles  well  below  the  static  stall  angle  (SSA),  in  the  neighborhood  of  and 

well  exceeding  the  time-averaged  dynamic  stall  angle  (DSA).  We  anticipated 

that  the  hysteresis  loops  would  be  different  from  those  of  an  oscillating 

12  13 

airfoil  as  derived  from  wind  tunnel  experiments.  *  Such  information  would 
shed  light  on  the  formation  of  the  hysteresis  in  the  three  different  flow 
regimes.  Comparison  will  be  made  with  laboratory  results  obtained  by  an 
oscillating  airfoil  in  a  steady  air  stream.  Finally,  we  will  present  the 
frequency  spectra  of  the  wind  components,  forces  and  moment.  The  spectral 
decomposition  of  these  quantities  would  provide  additional  information 
concerning  the  unsteady  aerodynamics  which  lead  to  lift  overshoot,  dynamic 
stall  and  potentially  to  stall  flutter. 


Time  Series  of  Ambient  and  Aerodynamic  Parameters 

Figures  2  through  4  show  the  time  series  of  the  three  wind  components,  U, 
V,  and  W,  the  lift,  drag,  pitching  moment  and  angles  of  attack  for  one  of  the 
runs  with  a  *  3.4°.  The  W  component  is  shifted  by  -6  m/ s  to  improve  the 
readability  of  the  figures.  The  Reynolds  number  and  the  turbulence  intensity 


are  Rc  ■  318,000  and  0.12,  respectively.  The  averaged  angle  of  attack  a 
corresponds  to  the  condition  at  the  minimum  drag  (D  =*  3.65  N  and  Cp  =*  0.025) 
[see  Table  1  and  Reference  8],  In  general,  the  lift,  drag  and  pitching  moment 
depend  on  both  U  and  a,  whereas  the  force  coefficients  depend  only  on  a  within 
a  narrow  range  of  Rc>  At  this  small  a,  increases  linearly  with  a  while  Cp 
has  the  minimum  value.  The  lift  and  drag  increase  whereas  the  pitching  moment 


decreases  with  increasing  U  and  a.  From  the  figures,  the  lift  and  pitching 
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Table  1.  Time-Averaged  Ambient  and  Aerodynamic  Parameters  for  the  High-R  Series  (Cont 
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moment  correlate  reasonably  well  with  both  U  and  a,  whereas  the  drag  correlates 
better  with  a.  It  is  evident  from  the  figures  that  the  lift  and  the  pitching 
moment  follow  the  variations  of  the  U  component  for  both  high-  and  low-frequency 
fluctuations.  The  W  component  or  a  and  the  drag,  on  the  other  hand,  consist 
of  only  relatively  high  frequency  fluctuations.  From  now  on,  we  emphasize  only 
the  dependence  of  the  aerodynamic  parameters  on  a,  which  governs  the  unsteady 
aerodynamics  and  thus  the  dynamic  stall  phenomenon. 

Note  that  the  fluctuations  in  the  drag  result  in  instantaneous  negative 
drag  (i.e.,  propulsion).  The  drag  becomes  negative  only  when  the  vertical  wind 
component  or  the  angle  of  attack  become  negative,  although  the  time-averaged 
drag  is  positive.  The  exact  mechanism  that  generates  the  instantaneous  pro¬ 
pulsion  is  not  known.  Under  certain  circumstances,  however,  the  unsteady  wind 
field  could  create  a  situation  similar  to  that  of  a  fish  which  manages  its  tail 
and  fins  to  propel  itself  through  unsteady  flows  around  its  body. 

Figures  5  through  7  show  the  results  for  a  “  19.9°,  u'/u  “  0.156,  and 
Rc  =  420,900.  For  such  large  a,  the  wing  is  in  and  out  of  dynamic  stall 
repeatedly  while  the  aerodynamic  parameters  vary  drastically,  as  demonstrated 
by  the  highly  fluctuating  time  series  that  consist  of  many  sharp  spikes.  From 

Q 

Liu,  the  time-averaged  DSA  for  the  Wortmann  wing  is  estimated  to  be  roughly 
19°  whereas  the  SSA  of  the  Wortmann  wing  is  between  14°  to  15° ,  as  derived 
from  wind  tunnel  results . ^ ^  Consider  the  part  of  a  typical  cycle  when  the 
W  component,  and  thus  a,  has  reached  the  minimum  and  begins  to  increase  again. 

At  this  point,  the  angle  of  attack  is  below  the  time-averaged  DSA  or  even  the 
SSA.  The  lift  and  drag  increase  with  increasing  a,  whereas  the  pitching  moment 
changes  insignificantly  with  a.  The  same  trend  continues  when  a  exceeds  the 
SSA  but  is  below  the  DSA. 

Drastic  changes  in  the  lift,  drag  and  pitching  moment  take  place  when  a 
continues  to  Increase  and  exceeds  the  DSA.  The  lift  experiences  a  sudden  drop 
as  a  exceeds  the  DSA.  When  the  amplitude  of  the  fluctuation  in  a  is  large, 
the  sudden  drop  in  the  lift  results  in  generating  higher  harmonics  in  the  time 
series.  The  lift  remains  low  for  the  rest  of  the  cycle.  Recovery  of  the  lift 
does  not  take  place  until  the  end  of  the  cycle  as  a  has  gone  through  the  maxi¬ 
mum  and  returned  to  the  initial  value  below  the  DSA.  The  drag  time  series  shows 
different  behaviors  from  that  of  the  lift.  In  comparison,  the  drag  increases 
with  increasing  a;  it  continues  to  increase  when  a  exceeds  the  DSA.  The  large 
spikes  in  the  time  series  of  the  drag  and  a  are  in  phase,  but  those  of  the  lift 
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and  drag  are  not  in  phase.  The  large  spikes  in  lift  are  accompanied  by  the 
emergence  of  higher  harmonics  in  the  time  series.  The  reduction  of  the  drag  to 
its  original  value  takes  place  at  the  end  of  the  cycle  when  a  returns  to  its 
initial  value  below  the  DSA.  There  also  exists  a  hysteresis  loop  for  the  drag. 

The  pitching  moment  does  not  change  considerably  for  a  between  3.4°  and  the 
DSA,  as  can  be  seen  from  the  near  constant  in  that  range  of  a  (Figure  8 
in  Reference  8).  As  the  DSA  is  exceeded,  the  pitching  moment  drops  rapidly,  in 
the  form  of  negative  spikes  shown  in  the  time  series.  Note  that  the 
contribution  of  the  lift  is  larger  than  that  of  the  drag,  for  the  Wortmann  wing, 
due  to  a  relatively  large  moment  arm  of  the  former  (to  the  quarter  chord). 

Figures  8  through  10  illustrate  the  results  for  a  *  23.0°,  u'/U  “  0.13,  and 

Rc  °  360,000.  For  this  run,  a  is  about  3°  higher  than  the  time-averaged  DSA. 

Comparison  of  the  aerodynamic  parameters  shows  that  there  are  more  large- 

amplitude  jumps  in  the  lift  time  series  for  a  =*  23.0°  than  for  a  ■  19.9°.  The 

corresponding  average  and  root-mean-square  (rms)  lift  coefficients  are 

CT  -  1.72  and  CT '  =  0.31  for  a  **  23.0°  and  CT  -  1.86  and  CT  '  “  0.23  for 
■Li  Jj  Li  Li 

a  ■  19.9°  (see  Table  2).  The  drag  and  pitching  moment  time  series  of  the  two  runs 
show  similar  characteristics.  From  Table  2,  the  average  and  rms  coefficients 
are  CD  “  0.47,  CD*  ■  0.21,  ■  -0.27  and  C^'  ■  0.066  for  a  •  23.0°  and 

CQ  =  0.31,  Op'  ■  0.16,  =  -0.23  and  C^’  =  0.050  for  0  ■  19.9°.  The  above 

comparison  indicates  that  the  wing  is,  on  the  average,  dynamically  stalled 
(light  stall,  see  Reference  10)  more  often  for  a  ■  23.0°  than  for  a  ■  19.9°, 
as  anticipated. 

Figures  11  through  13  show  the  results  for  a  ■  29.3°,  u'/U  =  0.12,  and 
R£  =*  331,000.  The  Wortmann  wing  is  in  deep  stall  (see  Reference  20)  most  of 
the  time.  The  time  series  of  the  aerodynamic  parameters  show  a  significant 
reduction  in  the  number  density  of  large-amplitude  spikes.  The  coefficients 
are  CL  -  1.48,  C  •  -  0.24,  CQ  -  0.76,  C  '  “  0.19,  -  -0.32  and  -  0.052, 

respectively.  Here,  the  lift  coefficient  reduces  significantly,  whereas  the 
drag  coefficient  continues  to  increase,  from  their  values  before  dynamic  stall 
takes  place.  Some  reductions  in  the  rms  coefficients  are  evident,  but  the 
values  of  these  coefficients  remain  reasonably  high. 

In  summary,  the  results  presented  here,  together  with  those  presented  pre- 
8 

viously  by  Liu,  show  that  the  unsteady  wind  field  at  6  u  above  the  ground 
reduces  the  minimum  drag  coefficient  of  the  Wortmann  wing  (from  about  0.05  to 
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Series  of  Wind  Components  (0.1-s  averages)  Measured  at  the  Ellensburg  Airport 
3.0*,  Rc  -  360,000. 
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Figure  10.  Tine  Series  of  Lift,  Drag  and  Pitching  Moaent  (0.1-s  averages) 
a  -  23.0°,  Rc  -  360,000. 
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0.025)  at  small  a  and  Increases  the  maximum  lift  coefficient  (from  about  1.6  to 

1.9)  at  large  a  beyond  the  SSA  and  below  the  DSA.  Wind  tunnel  data  (see,  for 

example.  Reference  4)  using  the  same  wing  model  in  a  steady  air  flow  were  used 

as  a  reference  for  the  present  full-scale  results.  The  SSA  and  the  DSA  are 

17  21  8 

measured  to  be  about  14°  to  15°  *  and  19°,  respectively.  For  a  exceeding 

the  dynamic  stall  angle,  referred  to  as  light  stall, ^  the  time  series  of  the 
aerodynamic  parameters  are  characterized  by  highly  fluctuating  signals  with 
large-amplitude  spikes.  The  average  CT  decreases  slightly,  and  C  continues  to 
increase  from  the  values  before  dynamic  stall.  C,  remains  nearly  constant  from 
a  ■  0°  to  15°  and  begins  to  decrease  for  a  exceeding  the  DSA.  The  rms  coeffi¬ 
cients  of  the  fluctuating  force  are  usually  highest  under  light  stall  condi¬ 
tions.  For  deep  stall  with  a  greatly  exceeding  the  DSA,  CT  and  C,.  continue  to 
decrease,  whereas  CQ  increases.  The  rms  coefficients  decrease  somewhat  but 
their  values  remain  relatively  high. 
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Hysteresis  Loops  of  Aerodynamic  Coefficients 

The  unsteady  coefficients  of  lift,  drag  and  pitching  moment  of  nine 
selected  cycles  from  the  time  series  with  relatively  large  amplitudes  of  the 
fluctuating  a  are  shown  in  Figures  14  through  22.  In  each  group  of  figures, 
we  Illustrate  the  lift,  drag  and  moment  coefficients  as  the  ordinate  versus 
the  Instantaneous  angle  of  attack  as  the  abscissa  for  three  individual  cycles 
with  a  =  3.4°,  19.9°  and  29.3°.  The  corresponding  time-averaged  coefficients 
are  plotted  as  a  dashed  curve  in  each  figure.  The  shapes  of  the  loops  appear 
to  be  considerably  different  for  each  Individual  cycle.  However,  there  are 
some  common,  organized  features  in  these  cycles.  Efforts  are  made  to  examine 
the  hysteresis  loops  in  these  coefficients  and  to  compare  then  with  wind  tunnel 
results  for  an  oscillating  airfoil.  It  is  important  to  bear  in  mind  that  there 
are  basic  differences  in  the  wind  tunnel  and  EATS  experiments.  The  unsteadi¬ 
ness  of  the  present  experiments  is  due  entirely  to  the  wind  field  consisting  of 

eddies  with  a  wide  spectrum  of  sizes  and  frequencies.  The  wind  tunnel  experi- 
12  13 

ments,  *  on  the  other  hand,  were  conducted  with  a  pitching  or  plunging 
airfoil  oscillating  at  a  single  frequency  in  a  steady  air  flow  with  a  certain 
turbulence  level  inherent  to  the  wind  tunnels.  Specifically,  the  unsteadiness 
of  the  present  experiments  is  due  to  the  variations  in  all  three  wind  compo¬ 
nents.  The  net  result  is  to  introduce  a  combined  unsteady  motion  that  includes 
streamwise  translating,  spanwise  sweeping,  pitching  and  plunging  components 
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simultaneously  and  that  is  three-dimensional  in  nature.  All  these  unsteady 
components  influence  the  performance  of  an  airfoil  (see  Reference  8).  In 
addition,  the  Reynolds  number  varies  within  a  relatively  narrow  range. 

Although  complicated,  these  are  the  realistic  conditions  to  be  encountered  by 
a  low-Reynolds-number,  low-altitude  RPV  deployed  near  the  ocean  surface. 
Consequently,  we  anticipate  that  the  shape  of  the  hysteresis  loop  will  vary 
significantly  in  the  unsteady  but  turbulent  wind  field  in  comparison  with 
those  of  the  single-frequency  oscillating  airfoil. 

Figures  14,  17  and  20  show  three  sets  of  hysteresis  loops  of  C  for 

1j 

a  “  3.4°,  19.9°  and  29.3°.  To  differentiate  the  adjacent  loops,  the  one  for 
a  =  19.9°  is  coded  with  dots.  Small  arrows  are  added  to  direct  the  time  paths 
of  the  loops.  The  loops  always  start  from  the  end  with  increasing  a  and  end 
with  decreasing  a,  as  directed  by  the  small  arrows.  The  shapes  of  the  loops 
appear  to  be  vastly  different  for  individual  cycles,  but  they  do  exhibit  some 
common,  organized  features.  The  loops  for  C  at  small  a  are  relatively  narrow, 

Li 

similar  to  those  of  an  oscillating  airfoil.  The  shape  is,  however,  more  com¬ 
plicated,  with  the  coexistence  of  miniloops,  and  it  varies  with  different  cycles 
For  a  >  19.9°,  the  instantaneous  CL  continues  to  increase  after  a  passes  the 
time-averaged  DSA  (i.e.,  19°),  which  is  determined  from  the  time-averaged 
versus  a  curve  (Figure  6  in  Reference  8).  In  some  cases,  the  sudden  drop  in  CT 

Li 

occurs  at  a  value  of  a  way  beyond  25°.  Evidently,  there  is  a  time  lag  between 
the  change  in  the  wind  field  and  the  flow  field  around  the  wing  as  It  is  re¬ 
sponding  and  adjusting  to  the  unsteady  wind  field.  This  lag  has  been  observed 

12  13 

and  documented  in  laboratory  results  of  an  oscillating  airfoil.  *  Once  in  a 
deep  stall,  CL  remains  low  as  a  reaches  its  maximum  value  and  begins  to  decrease. 
For  the  loop  with  a  “  29.3°,  the  wing  Is  in  deep  stall  most  of  the  time  and 
may  remain  low  throughout  the  entire  cycle.  The  lift  is  usually  higher  during 
the  first  half  of  the  cycle  with  increasing  a  than  it  is  during  the  second  half 
with  decreasing  a.  The  lift  often  restores  to  the  neighborhood  of  the  initial 
value  near  the  end  of  the  cycle  as  the  minimum  o  is  approached,  with  some 
exceptions.  The  loop  for  a  ■  29.3°  is  very  disorganized,  which  reflects  the 
highly  separated  flow  on  the  upper  surface  of  the  wing  in  deep  stall  regime. 

Figures  15,  18  and  20  show  three  sets  of  hysteresis  loops  of  for 
a  M  3.4°,  19.9°  and  29.3°.  Many  miniloops  are  present  within  the  major  loops 
of  each  Individual  cycle.  Again,  the  loops  for  the  small  a  are  relatively 
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narrow,  as  expected.  The  loops  for  the  large  angles  not  only  have  large 
amplitudes  but  also  show  complex  patterns.  Each  loop  deviates  significantly 
from  the  time-averaged  curve  (dashed).  For  a  *  19.9°,  there  is  a  sharp  rise 
in  Cp  during  the  increasing  a  half  of  the  cycle  for  all  three  of  the  examples. 
Such  a  sharp  rise  corresponds  to  either  an  abrupt  increase  in  CL  or  to  the 
shifting  into  stall  onset  and  light  dynamic  stall  regimes  as  a  approaches  and 
passes  through  the  DSA.  Note  that  the  bulk  of  is  composed  of  two  components, 
the  coefficients  of  the  form  drag  and  the  induced  drag;  the  latter  is  propor¬ 
tional  to  C  .  In  the  light  stall  regime,  C  levels  off  and  begins  to  drop  off, 

L  Li 

which  reduces  the  contribution  of  the  induced  drag  to  while  the  form  drag 
continues  to  increase.  In  the  deep  stall  regime,  drops  off  rapidly  and  re¬ 
mains  low,  further  reducing  the  contribution  of  the  induced  drag.  The  form 
drag  increases  significantly  as  a  result  of  the  formation  of  a  large  separated 
zone  on  the  upper  wing  surface,  which  maintains  the  increase  in  but  with  a 
more  gentle  slope  due  to  reduced  contribution  from  the  induced  drag. 

The  hysteresis  loops  for  are  illustrated  in  Figures  16,  19,  and  22. 

Note  that  is  negative  for  a  nose-down  orientation.  Knowing  that  each  loop 
deviates  from  others  significantly  due  to  the  random  nature  of  the  cycle,  we 
attempt  in  the  following  to  describe  a  general  trend  of  the  loops.  Because 
both  CT  and  C_  contribute  to  C..,  the  steepest  drop  or  rise  of  C,  corresponds 
to  the  situation  when  and  CD  rise  or  fall  simultaneously.  Therefore, 
drops  abruptly  when  a  approaches  and  passes  through  the  DSA  before  the  stall 
onset,  while  it  rises  rapidly  when  a  has  passed  the  DSA  and  the  wing  is  experi¬ 
encing  the  deep  dynamic  stall  process.  drops  again  when  the  wing  remains 
stalled  dynamically  during  which  is  low  and  rises  due  to  the  Increase  in 
the  form  drag.  Finally,  rises  and  returns  to  roughly  the  initial  value  at 
the  end  of  the  cycle  as  a  passes  the  maximum  angle  during  that  cycle  and  begins 
to  decrease  again.  For  a  cycle  with  an  amplitude  of  5°  to  6°,  the  excursion  of 
C,^  could  be  0.2  or  larger  for  a  “  19.9°  and  29.3°. 

Of  all  the  C,^  loops,  only  the  one  for  a  ■  19.9°  shows  a  negative  damping 
loop  (clockwise).  At  that  angle,  the  wing  is  in  and  out  of  dynamic  stall,  a 
favorable  condition  for  negative  damping.  During  that  cycle,  the  area  of  the 
clockwise  loop  is  considerably  larger  than  that  of  the  counterclockwise  loop. 
Under  this  condition,  the  wing  could  experience  stall  flutter  provided  its  struc¬ 
tural  natural  frequency  is  close  to  the  dominant  frequency  of  the  wind  field. 
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The  frequency  spectra  of  U,  W,  a,  L,  D  and  M  for  four  runs  with  a  “  3.4° 
(rectangles),  19.9°  (pluses),  23.0°  (diamonds)  and  29.3°  (triangles)  are  shown 
in  Figures  23  through  28.  The  spectra  were  estimated  by  calculating  the  power 
spectra  of  the  individual  time  series  averaged  over  0.1  second.  There  were  a 
total  of  900  data  points.  First,  we  removed  the  linear  trend  from  the  time 
series  and  then  added  62  zeros  to  each  end  of  the  time  series  to  make  up  a 
total  of  1024  points  for  evaluating  the  power  spectra.  A  cosine  window  as 
described  by  Bendat  and  Piersol  was  used  to  reduce  leakage  outside  the  fre¬ 
quency  band  f  ■  +  1/T,  where  T  is  the  total  time  or  102.3  seconds.  The  power 
spectra  were  evaluated  by  using  the  algorithm  available  in  the  ASYSTANT 
scientific  software  (Macmillan  Software  Company)  for  IBM  PC/AT  or  compatibles. 

To  reduce  the  spread  of  the  spectral  density,  a  running  mean  algorithm  was 
applied  by  averaging  one  to  nine  points  of  the  frequency  distribution,  with 
gradual  increase  in  the  number  of  points  from  the  low-  to  high-frequency  end. 

The  spectra  of  the  U  and  W  wind  components  for  the  four  runs  (Figures  23 
and  24)  show  that  the  shapes  of  the  spectra  are  quite  similar.  The  difference 
in  the  spectral  density  is  mainly  due  to  the  variations  in  the  turbulence 
intensity  of  the  individual  runs,  as  anticipated.  From  Table  1,  we  determined 
that  the  turbulence  intensities  u'/U  for  the  four  runs  are  0.12,  0.16,  0.13  and 
0.12,  respectively.  Figures  23  and  24  show  that  the  amplitude  of  the  spectral 
density  is,  on  the  average,  proportional  to  u'/U.  The  same  characteristics 
are  observed  for  the  spectra  of  a  (Figure  25),  which  are,  of  course,  closely 
related  to  those  of  U  and  W.  The  U  spectra  show  a  -5/3  slope  for  about  one 
decade  up  to  f  a  1  to  2  Hz.  Beyond  2  Hz,  the  spectral  density  rolls  off  at  a 
slope  slightly  steeper  than  -5/3.  As  discussed  in  Liu,  the  steep  rolloff  is 
due  to  an  Insufficient  frequency  response  of  the  Gill  propeller  anemometer 
whose  -3  dB  points  is  about  2  Hz  at  U  =*  8  m/s .  The  W  spectra  (Figure  24)  show 
a  slightly  narrower  -5/3  region  than  the  U  spectra.  Beyond  2  Hz,  the  rolloff 
is  slightly  steeper  for  the  W  spectra  than  for  the  U  spectra.  Again,  as 

3 

discussed  by  Liu,  the  response  in  the  W  component  of  the  Gill  propeller  is 
worse  due  to  the  small  value  of  W.  A  partial  remedy  was  made  by  tilting  the 
anemometer  at  an  angle  of  15°  to  Improve  the  frequency  response  of  W.  The  a 
spectra  show  a  shape  similar  to  that  of  the  W  spectra,  because  a  ■  arctan  (W/U). 
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The  L,  D  and  M  spectra  (Figures  26  through  28)  all  display  a  region  of 

-5/3  slope  up  to  4  Hz,  which  is  higher  than  those  of  the  U  and  W  spectra.  As 
0 

described  by  Liu,  the  frequency  response  of  the  aerodynamic  balance  was 
significantly  higher  than  that  of  the  Gill  anemometers.  To  match  the  frequency 
response  of  the  two  instruments,  the  outputs  of  the  aerodynamic  balance  were 
low-pass  filtered  at  4  to  5  Hz  (-3  dB).  This  matching  was  essential  to  avoid 
artificial  spikes  in  the  spectra  of  C. ,  C_  and  C.  due  to  the  above  mismatch, 
as  the  coefficients  were  derived  from  simultaneous  measurements  of  the  wind  and 
force  components.  It  should  be  pointed  out  that  the  L,  D  and  M  spectra  depend 
on  the  variations  in  U  and  W  or  a.  Attempts  will  be  made  to  isolate  the 
dependence  of  these  parameters  by  correcting  the  spectra  with  the  turbulence 
intensity.  Without  the  correction,  the  effects  of  dynamic  stall  on  the  spectra 
cannot  be  assessed  due  to  the  interdependence  of  the  spectra  on  U  and  a.  In 
particular,  the  L  spectra  (Figure  26)  tend  to  be  a  stronger  function  of  U  than 
of  a.  The  effects  of  the  dynamic  stall  on  the  L  spectra  are  therefore  masked 
by  those  of  U.  The  D  and  M  spectra  (Figures  27  and  28),  on  the  other  hand, 
show  a  stronger  dependence  on  a  than  on  U.  For  example,  the  spectral  densities 
for  a  =  3.4°  (rectangles)  are  significantly  lower  than  those  for  a  ■  19.9°, 
23.0°  and  29.3°.  For  a  >  19°,  however,  the  effects  of  the  dynamic  stall  again 
cannot  be  isolated  from  those  of  U.  Note  that  there  are  two  distinct  spikes  at 
f  »  2.7  and  3.5  Hz  in  the  D  spectra  for  a  =*  3.4°.  The  spikes  persist,  but  the 
frequencies  shift  slightly  for  the  large  angles.  No  such  spikes  are  present  in 
either  the  L  or  M  spectra. 

To  isolate  the  dependence  of  a  from  that  of  U,  we  corrected  the  U,  W  and 
a  spectra  by  dividing  them  with  the  square  of  the  turbulence  intensity  or 

2  _ l 

u'  /U  .  Such  a  correcting  factor  was  selected  based  on  the  physically 

2 

consistent  relation  /^(f)  df  *  u'  where  ^  is  the  spectral  density  of  U. 
Figures  29  through  31  illustrate  the  corrected  U,  W  and  a  spectra.  It  is 
evident  that  the  correction  has  collapsed  the  spectra  and  has  reduced  the  de¬ 
pendence  of  U  significantly.  Similarly,  we  corrected  the  L,  D  and  M  spectral 

densities  by  dividing  them  with  the  u'  /U  ,  which  is  again  physically  con- 

2 

sistent  because  L  <*  U  .  Figures  32  through  34  illustrate  the  corrected  L, 

D  and  M  spectra. 
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From  the  corrected  L  spectra  (Figure  31),  we  observe  that  the  spectral 
densities  for  a  3  3.4°  is  generally  lower  at  low  (f  <  0.3  Hz)  and  high 
(f  >  0.3  Hz)  frequencies  than  those  for  a  3  23.0°  and  29.3°.  They  become 
roughly  the  same  for  medium  frequencies  (0.3  Hz  <  f  <1.3  Hz).  This  reflects 
the  presence  of  low-  and  high-frequency  fluctuations  beyond  the  stall  onset 
angles  for  L,  D  and  M.  Note  that  the  runs  for  a  3  23.0°  and  29.3°  correspond 
to  the  wing  in  the  light  and  deep  dynamic  stall  regimes,  respectively,  most  of 
the  time.  For  a  3  19.9°,  with  the  wing  operating  constantly  in  and  out  of 
dynamic  stall  (stall  onset),  the  spectral  densities  are  the  highest  at  low 
frjq>.  .ucies  (f  <  0.4  Hz).  For  f  >  0.4  Hz,  the  spectral  densities  are  inbetween 
‘I' ’St  for  a  3  3.4°  and  23.0°.  The  above  observations  indicate  that  the  L 
spectra  are  generally  low  for  all  but  the  medium  frequencies.  In  the  vicinity 
of  stall  onset,  the  L  spectral  densities  increase  the  most  at  low  frequencies, 
whereas  light  and  deep  stall  promote  high-frequency  fluctuations  in  the 
presence  of  highly  separated  flows  on  the  upper  wing  surface. 

The  corrected  D  spectra  (Figure  33)  show  more  dramatic  differences.  For 

a  3  3.4°,  the  spectral  density  is  consistently  lower  than  those  of  the 

others  for  all  frequencies.  The  presence  of  two  distinct  spikes  at  f  *  2.7 

and  3.6  Hz  is  evident,  but  their  source  is  not  known  at  present.  For 

a  3  19.9°,  the  spectral  density  is  significantly  higher  than  that  for 

a  3  3.4°  but  consistently  lower  than  those  for  a  3  23.0°  and  29.3°.  The 

spectral  densities  for  a  3  23.0°  and  29.3°  are  quite  similar;  the  latter  is 

higher  with  more  distinct  spikes  at  high  frequencies.  This  Is  consistent  with 

the  finding  that  the  hysteresis  loops  have  relatively  large  excursions  as  the 

wing  is  operating  in  and  out  of  the  deep  dynamic  stall  regime.  Furthermore, 

we  anticipate  that  the  flow  for  a  3  23.0°  is  not  as  strongly  separated  as 

that  for  a  3  29.3°,  as  CT  drops  off  and  remains  low  whereas  C  continues  to 

■L*  1) 

rise  when  a  Increases  from  23.0°  to  29.3°. 

The  corrected  M  spectra  (Figure  34)  have  shapes  between  those  of  the  L  and 
D  spectra  because  both  L  and  D  contribute  to  form  M.  One  of  the  most  distinct 
differences  in  the  M  spectra  is  that,  for  a  3  3.4°,  they  have  a  slope  con¬ 
siderably  steeper  than  those  of  the  rest  (roughly  -5/3  slope).  The  two 
distinct  spikes  at  f  »  2.7  and  3.6  Hz,  as  observed  in  the  D  spectra,  are  no 
longer  present  in  the  M  spectra.  At  the  high-frequency  end,  the  spectra  for 
a  3  29.3°  have  the  highest  spectral  density,  followed  sequentially  by  those 
for  a  -  23.0°,  19. 9 1  and  3.4°. 
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Comparison  of  the  L,  D  and  M  spectra  indicates  that  the  unsteady  wind 
field  has  the  most  impact  on  D  followed  by  M  and  has  the  least  impact  on  L,  as 
far  as  the  high-frequency  components  are  concerned.  The  dynamic  stall  has 
significant  effects  on  these  spectra  at  frequency  greater  than  about  2  Hz 
where  noticeable  increase  in  the  spectral  density  is  observed.  When  the  wing 
is  going  through  light  and  deep  dynamic  stall,  especially  the  latter,  the 
spectra  becomes  spiky  with  relatively  well-defined  dominant  frequencies.  The 
spiky  spectra  are  a  signature  of  the  presence  of  large  hysteresis  loops. 
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4.  SUMMARY 

An  environmental  aerodynamic  test  system  (EATS)  was  designed  and  assembled 
to  investigate  the  effects  of  atmospheric  unsteadiness  (turbulence  and  gust)  on 
the  performance  of  a  full-scale  Wortmann  FX-63-137  wing  (smoothed  by  Eppler) 
(see  Section  2  and  also  References  8  and  16).  Several  series  of  experiments 
were  conducted  under  a  "park-mode"  configuration  by  directing  the  wing  into 
the  prevailing  wind  at  windy  sites  (see  Section  3  and  Reference  8).  The 
time-averaged  aerodynamic  parameters  derived  from  the  full-scale  experiments 
differ  considerably  from  those  of  the  wind  tunnel  counterparts  under  steady 
conditions.  The  unsteady  wind  field  was  demonstrated  to  have  significantly 

g 

favorable  effects  on  the  performance  of  the  Wortmann  wing.  The  most  important 
unsteady  effects  pertain  to  the  lift  overshoot,  stall  delay,  reduction  of  drag 
and  endurance  enhancement.  The  lift  overshoot  and  stall  delay  are  consistent 
with  wind  tunnel  observations  for  an  oscillating  airfoil.^ 

To  understand  the  above  unsteady  effects,  further  analysis  of  the 
full-scale  data  was  conducted  with  emphasis  on  examining  the  time  series,  the 
unsteady  loops  and  the  frequency  spectra  of  the  ambient  and  aerodynamic 
parameters.  The  salient  points  derived  from  the  analysis  are  summarized  below. 

(1)  For  small  a  =<  3.4°,  which  is  in  the  minimum  drag  regime,  the  time 
series  of  D  shows  instantaneous  negative  values  (i.e.,  propulsion), 

g 

although  the  time-averaged  drag  is  positive.  The  exact  mechanism 
that  generates  the  instantaneous  propulsion  is  not  known. 
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(2)  For  large  a  in  the  neighborhood  of  the  dynamic  stall  angle,  the  time 
series  of  L,  D  and  M  contain  some  spikes  with  large  fluctuations.  The 
unsteady  wind  field  changes  the  instantaneous  angle  of  attack  of  the 
Wortmann  wing  and  causes  it  to  experience  the  dynamic  stall 
phenomenon,  which  is  drastically  different  from  the  static  stall 
phenomenon.  During  any  one  of  the  random  cycles  with  a  sufficiently 
large-amplitude  fluctuation  in  the  W  component,  the  wing,  which  is 
initially  not  stalled,  may  go  through  all  the  dynamic  stall  regimes 
(stall  onset,  light  and  deep  stall): 

o  Before  stall  onset,  significant  lift  overshoot  was  observed, 

resulting  in  a  large  increase  in  C  (from  1.6  to  over  2.0).  The 

•Li 

instantaneous  stall  onset  takes  place  after  a  greatly  exceeds 


8 

the  time-averaged  dynamic  stall  angle  (about  19°).  This  stall 
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delay  is  due  to  the  time  lag  for  the  adjustment  of  the  unsteady 
flow  field  around  the  wing  in  response  to  the  change  of  the  wind 
field  (Figures  14  through  22).  The  lift  overshoot  leads  to  a 
considerable  increase  in  the  drag  due  to  the  increase  in  both  the 
form  and  induced  drag.  The  pitching  moment  begins  to  drop  off 
sharply  as  the  lift  overshoots  its  maximum  static  value. 

At  stall  onset,  the  lift  levels  off  and  drops  slightly.  The 
pitching  moment  shows  similar  behavior.  The  drag  appears  to  rise 
continuously  due  to  the  increase  in  the  form  drag,  although  the 
contribution  of  the  induced  drag  decreases. 

Further  increase  in  a  causes  the  wing  to  enter  the  light  stall 
regime  with  the  lift  dropping  off,  the  pitching  moment  rising  and 
the  drag  dropping  off  momentarily. 

Negative  aerodynamic  damping  is  experienced  by  the  wing  as  it  is 
in  and  out  of  stall.  Stall  flutter  may  be  initiated  provided  the 
natural  structural  frequency  of  the  wing  is  close  to  the  dominant 
frequency  of  the  wind  field. 

As  the  wing  enters  the  deep  stall  regime,  the  lift  drops  below  the 
maximum  time-averaged  value.  During  that  time,  the  thickness  of 
the  separation  region  increases  to  the  scale  of  the  chord.  As  a 
result,  the  form  drag  rises  sharply.  The  total  drag  rises  more 
moderately  due  to  the  decrease  in  the  induced  drag.  Meanwhile, 
the  pitching  moment  drops  to  its  minimum  value. 

As  a  reaches  its  maximum  and  begins  to  decrease,  the  lift,  drag 
and  pitching  moment  begin  their  recovery.  At  the  end  of  the  cycle, 
they  restore  to  roughly  the  initial  values.  The  paths  during  the 
two  halves  of  the  cycle  differ  significantly,  which  forms  the 
hysteresis  loops.  Due  to  the  randomness  of  the  wind  field,  the 
hysteresis  loops  consist  of  many  miniloops  and  are  much  more 
complicated  than  those  of  an  oscillating  airfoil  with  a  single 
frequency.  The  dynamic  stall  process  is  presented  in  terms  of  the 
response  in  the  forces  and  moment  as  can  be  seen  in  both  the  time 
series  and  the  unsteady  loops  (Figures  2  through  22.) 


(3)  The  frequency  spectra  of  the  ambient  parameters  demonstrate  that  the 

spectral  shapes  are  similar  for  the  four  runs  selected  for  detailed 

investigation,  i.e.,  with  "a  ■  3.4°,  19.9°,  23.0°  and  29.3°.  The 

spectral  densities  for  U,  W  and  a  are  proportional  to  U  and  a. 

The  spectra  collapse  reasonably  well  when  they  are  corrected  for  by 

2  2 

dividing  the  square  of  the  turbulence  intensity  or  u'  /U  .  Such 

physically  consistent  correction  effectively  removes  the  dependence  of 

the  spectra  on  U.  Similarly,  the  spectra  of  the  aerodynamic  parameters 

♦  _ % 

are  corrected  for  by  dividing  u’  /U  to  remove  their  dependence  on  U. 
Comparsion  of  the  corrected  spectra  of  the  aerodynamic  parameters  has 
shown  several  important  features: 

o  The  unsteady  wind  field  has  the  most  impact  on  the  drag  spectra 
followed  in  order  by  the  pitching  moment  and  lift  spectra. 

For  example,  the  spectral  density  of  the  drag  for  a  3.4°  is 
consistently  1  to  2.5  orders  lower  than  that  for  a  =  23.0°  or 


29.3°,  except  at  the  low-frequency  end  where  the  resolution  of  the 
spectral  analysis  is  less  adequate.  The  spectral  density  of  the 
pitching  moment  is  over  1  order  of  magnitude  lower  for  a  =  3.4° 
than  for  a  =*  23.0°  or  29.3°  for  frequencies  higher  than  1  Hz. 

The  shapes  of  the  L  and  M  spectra  show  a  steeper  slope  at  the 
high-frequency  end  for  a  =  3.4°  than  for  its  high  a  counter¬ 
parts,  especially  for  the  L  spectra, 
o  As  a  result  of  dynamic  stall,  the  spectral  densities  of  L,  D  and  M 
at  high  frequencies  increase  considerably  together  with  the 
introduction  of  large-amplitude  spikes  with  well-defined  dominant 
frequencies.  The  spikes,  especially  in  the  D  spectra,  are  most 
distinctive  for  a  ■  29.3°  when  the  wing  is  in  the  deep  stall 
regime  most  of  the  time.  The  distinct  spikes  are  the  signatures 
of  the  presence  of  large  hysteresis  loops. 
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